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Abstract 

An absolute luminosity can be measured in an accelerator by sweeping beams transversely 
across each other in the so called van der Meer scan. We prove that the method can be 
applied in the general case of arbitrary beam directions and a separation scan plane. A simple 
method to develop an image of the beam in its transverse plane from spatial distributions 
of interaction vertexes is also proposed. From the beam images one can determine their 
overlap and the absolute luminosity. This provides an alternative way of the luminosity 
measurement during van der Meer scan. 
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1. Van der Meer method for arbitrary beam velocities 

The luminosity of an accelerator is given by 

L = fNiN2K ■ J pi^^(f- Af, t)p2^^(f, t) d^fdt, (1) 



where K = \J {vi — V2Y — ^"^^2^^^ is the M0ller kinematic relativistic factor c is the 
speed of light, Ai^2 are the number of particles in the colliding bunches all moving with 
the common velocities ^1^2, / is the frequency of collisions and p^i\{r^t) are the normalized 
particle densities in the laboratory frame, so that j pf^{r,t) d^r = 1 at any time t. The 
absolute value of the luminosity or the cross section can be measured by separating the 
beams in the transverse plane by Af and by monitoring a collision rate as a function of 
Ar. This method was proposed by van der Meer more than 40 years ago and was originally 
proved in Q for arbitrary beam shapes and parallel beams tTi ||t'2. It was successfully applied 
with various modifications at ISR jH, [sj , RHIC j3] and recently at LHC jsf accelerators. It 
was often used in the approximation of Gaussian or double Gaussian beam shapes. At 
RHIC, for example, this allowed to take into account various corrections due to the so-called 
hourglass effect, the beam-beam deflection and the beam crossing angle. The latter alone, 
however, does not require any significant changes in the original van der Meer method. Since 
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we did not see any publication on this subject, in this section we present a proof of van der 
Meer formula in case of arbitrary beam crossing angle and beam shapes. It is applicable to 
the scans at LHC where hourglass and beam-beam effects are small [5j. 




Without loss of generality in Eq. ([T]) it is assumed that only the first beam is moved. We 
choose a coordinate system as shown in Fig. [1] with z axis along the direction Ai7 = V\—V2i 
X axis lying in the beam crossing plane and y axis perpendicular to x and z. Let's denote 
z- and x-components of the velocities as Vi^2z and v±, respectively, so that Vi^2 = v± + Vi^2z- 
The beam displacement plane is not necessarily perpendicular to z, and in general Ar has 
three components (Ax, Ay, Az). Its projection to x-y plane will be denoted by Arj_. For 
the particles uniformly moving with the velocities vi^2 the time evolution of their densities 
obeys the rule pf\{r, t) = p^i2{r — vi,2t, 0), therefore 
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where we changed the integration variables to r — v\t = (x — v\t^ y, z) 
z — Az — vizt = z' , z + V2zt = z", the corresponding Jacobians were 

_J{r^^ d{z,t) ^ 

d{r-v^t,t) ' d{z',z") ^' 

We also used the notation p^i2'^i^±) = / p'f'2(^J-5 0) dz for the particle density projections 
on the plane perpendicular to z \\ Av at t = 0. For the case Vi\\v2 when z || vi^2 we recover 
the usual formula 
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Following van der Meer method, Eq. ([2]) should be integrated over Ar. In three- 
dimensional space an equation of Ar* plane can be written in the form Ar* ■ n = A where 
A is some constant and n = {cos6x, cos By, cos6z) is the unit normal, O^^y^z are the angles 
between n and the axes. Since Eq. (|2]) does not depend on Az, an integration over Ar can 
be performed with the help of 5-function as 

dAx dAy d^Af^ 



d^Ar = 6{Af -n-A) dAx dAy dAz 
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In other words, the independence on Az allows to vary only Arj_ projection of Ar*, and 
1/ cos 6*2 appears as a difference between the area on Ar plane and its x-y projection. Inte- 
gration of the density product in Eq. ([2]) over dAx dAy = d'^Af± gives 
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and we obtain 
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For the process with the cross section a (including a reconstruction efficiency) the rate 
of events is given by 

R{Af)=a-L{Ar). (6) 

If /, A^^i 2 are measured during the scan, the monitoring of the rate R{Af) allows to determine 
the cross section and then the absolute luminosity 
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This is the generalized van der Meer formula which is valid for any crossing angle between 
the beams and for arbitrary displacement plane. Note, that if the latter is not perpendicular 
to z, cos ^2 7^ 1, Az component of the displacement affects the time of the interactions but 
not the luminosity. 

The term with the velocities Vi^2 in Eq- O is simply a 7-factor of the boost with the 
velocity v± (see Fig. [T]). Indeed, it can be calculated as follows 
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Its appearance here can be understood from the following alternative proof of Eq. ([7]). It is 
based on a relativistic invariance. We start from Eq. ([3]) used by van der Meer in his original 
paper 0] and valid for the coUinear beams V\ || V2- Integration over Arj_ gives 
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where "coU" superscript of Ar*/"'^ reminds us that the formula is vahd only in the frame 
where the beams are collinear. The element Ar'j^"" is invariant only under boosts along z 
which also preserve the condition v\ || v^- The ratio of two frequencies Rj j is an average 
number of interactions per collision, which does not depend on the choice of the coordinate 
system, o is also a relativistic invariant according to the cross section definition proposed 
by M0ller The only non-invariant quantity in Eq. ^ is the displacement Af^j^°^^. 

The general case v\ j]' (see Fig. [1]) can always be reduced to the collinear beams by 
making a boost with the velocity v^. It is easy to show that regardless of the laboratory 
z-components V\^2zi in the relativistically boosted frame the beams become parallel to z, 
so that Eq. ^ is valid. To transform it back to the laboratory frame we note that (i^Arj_ 
transforms in the same way as a transverse part dx dy of a relativistically invariant four- 
dimensional volume dx dy dz dt. Since dz is not affected by the boost while dt acquires a 
7-factor, we have d'^A'r± = (i^Arj^°"/7_L, so that the relativistically invariant generalization 
of Eq. is 

- = 7. d^Af. = 7. cos e. ! d^^r. (10) 
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Together with Eq. (|8]) this completes the proof of Eq. ([7]) . 

For the following discussion it is useful to introduce coordinates linked to the beams 
as shown in Fig. [T] They are denoted by "1,2" subscripts. zi^2 axes are chosen along the 
beams, Xi^2 axes lie in the beam crossing plane and other axes coincide, 2/1 = 1/2 = V- If ^1,2 
is the angle between zi^2 and z, we have 

a; = a;i_2 cos a;i^2 + -2i,2 sin ai^25 z = =F 3^1,2 sin ai^2 ± -21,2 cos ai^2 ■ (H) 

Let's assume that the distribution of particles in the transverse x-y plane is independent 
in X and y. In the beam's frame this means an independence in xi^2 cos 0:1,2 + -^1,2 sin ai,2 and 
y, which is usually ensured by an absence of xi_2 — y accelerator coupling and an indepen- 
dence of longitudinal (-^1,2) and transverse (y) particle distributions. With this assumption, 
the two-dimensional integral over Af±_ can be reduced to the product of one- dimensional in- 
tegrals along any Ax = Axq and Ay = Ay^ lines. Indeed, if the transverse densities factorize 
as 2 ("^5 y) = Pi,2 X (^) ■ Pi,'i'y (y)^ the luminosity and the rate according to Eqs. ([2]) and 
(Q can also be factorized into x- and ^/-dependent terms, R{Ax,Ay) = Rx{Ax) ■ Ry{Ay), 
and its integral over d'^Ar± can be expressed as 



R{Ax, Ay) dAx dAy = j R^^iAx) dAx j Ry{Ay) dAy 
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In Eq. f|T2l) we considered the simple case when the integration is performed over Ax- Ay 
plane, so that cos 6^2 = 1. After substituting Eq. f lT2|) into ffTOjl . the cross section can be 
written as 

jRjAx, Ayo) dAx ■ jRjAxo, Ay) dAy 

" fN,N2RiAxo,Ayo) ' ^ ' 
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The integrals in the enumerator of Eq. (|T2l) can be measured in two one-dimensional scans 
over Ax at fixed A?/o and vice versa. There is no need to make a full scan in Ax- Ay plane. 
Note, that to get larger rates it is advantageous to keep the beam separation (Axq, Ayo) 
small, but in general Eq. (1131) is valid for arbitrary {Axq, Ay^). 

In derivation of Eq. [T3] we do not use the fact that x axis lies in the crossing plane. If 
the transverse distributions p^i2'~^ are independent in two arbitrary directions x' and y' in 
Af± plane, the formula f|T3|) written in the primed coordinates still remains valid, and x' 
and y' may be chosen as scan axes. If they are not perpendicular but form an angle ax'y', 
one should also include the corresponding Jacobian sin a^/y'. Finally, the two scan axes 
may extend beyond the x-y plane if their projections still coincide with x', y' axes. Let's 
denote the displacements along such scan axes as Ax" , Ay" and their inclination angles to 
x-y plane as a^",?/"- Taking into account that dAx'iy') = cosax"{y")dAx"{y"), one obtains a 
more general formula 

jRjAx", Ay'i) dAx" ■ jRjAx'^, Ay") dAy" 
a = 7xsma.vcosa..cosa, f N,N,R{Axl Ay'^) ' ^^'^ 

2. Reconstruction of individual beam profiles 

The density of interaction vertexes accumulated during time AT is given by 

d^N^. 



dx dy dz 



fATNiN2Ka I Pi^{f- 6f- Af, t) P2^{f- 6f, t) V^''^{6r) dt d:'5r. (15) 



This expression is very similar to the luminosity formula ([T]) except it contains the convolu- 
tion with the experimental vertex resolution V^'^'^{5r) and there is no integration over d^f. 
In analogy with van der Meer method, we may integrate this equation over d'^Ar to drop 
out pif"^, and after deconvolution with V^''^^((5r) obtain the profile of the second beam p]f^ . 
A possible non-zero beam crossing angle complicates the formulas, but as it is shown below, 
it is always possible to reconstruct the profiles of both beams in their transverse planes 
without any simplifying assumptions on the beam shapes. 

It is convenient to express the particle density of the second beam in its own coordinate 
system (0:2, y, Z2) = r2 (see Fig. [1]) as 

P2(r2, t) = (f, t) = p2{4, Z2 - \V2\t, 0). (16) 

The transformation Eqs. (ITT]) define the relation between the densities pl^^ and p2 in the 
laboratory and in the second beam frames, respectively. In the last equation we distinguished 
transverse = {x2,y) and longitudinal Z2 coordinates and took into account that Z2 \\ V2- 
In addition, let's define the resolution viewed from the second beam frame as 

V2iSf2) = V'^\6f). (17) 

For the first beam we define R = f — Af — Vit, so that 

p\^^{f- 6f - Af, t) = p^f'i.R - 6f, 0). (18) 
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Since / p^f{R-5r, 0) d?R = 1, to drop out pl^'^ from Eq. ^ one needs a three-dimensional 
integration over R. The volume element d^R is invariant under rotations. The simplest is 
to write it in the coordinate system with x and y axes lying in the Ar* displacement plane 
and with z axis pointing along the unit normal n. In components, cPR = dX^ dY^ dZ^, 
where A superscript denotes coordinates in this system. Only in Eq. (|T5|) depends on 
Af, therefore regardless of other variables, integration over Ar is equivalent to integration 
of over dX^ dY^. The third coordinate = R ■ n is spanned when integrating over t 
due to \vi\t term in R. 

The density of the second beam, ^2(^2"? -^2 — |t^2|^, 0) also depends on t, however. To 
decouple pf"^ and p2, one may integrate in addition over Z2 of the reconstructed vertexes. 
Then it is possible to change the integration variables from dtdz2 to dZ^ dZ2, where Z2 = 
Z2 — \v2\t. dZ^ completes the integration over d^R, while the integration over dZ2 gives a 
projection transverse to the second beam 

P2{r^) = j P2(rt, Z2 - 6z2, 0) dZ2 = j P2(rt, z, 0) dz. (19) 

Since r in i? also depends on Z2 and dr/dz2 = V2/\v2\, the Jacobian of the variables substi- 
tution is 



d{t, Z2) 



d{Z^, Z2) 



\n ■ [vi - V2j 



^ (20) 



cos ^7 ■ \Av \ ' 



where cos 9z is the same angle as in Eq. ([5]) between n and z \\ {vi — €2)- 
Combining all pieces f lT5|) -f l20|) and ([8]) together, one has 

X j P2{4 - Z2 - 5z2, 0) dZ2 j V2{54, SZ2) d5z2 dHr^ = 

pii4 - s4)V2^i64) d'64. (21) 



Here were used the rotation invariance of the volume element dx dy dz = dx2 dy2 dz2 and the 
same for d^6f, and also defined 



^^(^f^) = / V2i5f^, 5z2) dSz2 = / V^^^{5r) dSz2, (22) 



which is the vertex resolution in the beam transverse plane. Integration of Eq. fl2T|) over 
f2 gives the total number of interactions and their rate consistent with Eq. ffTOl) . Unfolding 
Eq. fl2T]) with V^{6r2) gives the image of the second beam in its transverse plane P2"(^2')- 

From the very beginning it was assumed that only the first beam is moved in van der Meer 
scan. Instead, one can go to its "rest" frame where vice versa, only the second beam moves 



6 



by — Ar and the first is stable. By repeating the procedure above, one can then measure 
the transverse image of the first beam, as in this case the integration over Ar* and zi leads 
to the same Eq. (12T|) with the index 2 substituted by 1. It is interesting that two beam 
images can be obtained simultaneously from the same set of vertex distributions measured 
at various Ar points. They should be brought to the rest frame of the corresponding beam, 
i.e. aligned differently during summation, and projected to the beam transverse plane. 

If /, AT, or the reconstruction efficiency entering a change during the scan, they 
should be considered as functions of Ar* during the integration on the left side of Eq. fl2T|l . 
On practice, this means that the accumulated vertex distributions and their statistical errors 
should be reweighted according to the factor {fATNiN2<7)~^ at every scan point. 

One can determine not only transverse but the full three-dimensional images pi^2, if the 
detector is able to measure the time t of the interaction. Indeed, in this case Eq. (|T5|) may 
be integrated over Ar* and {f— vit) ■ n = Z^, while the integration over Z2 = Z2 — \v2\t 
may be omitted. This is possible since is directly measured. In addition to the vertex 
resolution V^'^^{6'r), we should include in the formulas the time resolution Vt{St), make a 
change t ^ t — 6t and integrate over 6t. As before, after the integration over (PArdZ^ or 
d^R the density p^i^ drops out and we obtain 

r\ f ^ ^ ^^YX 79a— *T/r— * — * 1 — *\ 

p\^^{f- 6f- V2{t - 6t), 0) V^^'^iSr) Vt{6t) d^fdSt. (23) 

The distribution of vertexes in (x2, 2/2, -^2 — W2\t) space, accumulated during the scan and 
deconvolved with the resolutions determine the density p2- Note, that the time resolution 
affects the imaging in the longitudinal Z2 \\ V2 direction. With an infinitely poor resolution 
in Z2 — \v2\t measurement, Eq. f l23|) becomes effectively equivalent to Eq. fl2T]) . 

Now let's consider the case when the particle distributions in y direction and in x-z 
crossing plane are independent. In the following, x,z- and ^/-dependent parts of a func- 
tion will be denoted by the corresponding subscripts, so that we have pY2i^y ^5 0) = 
Pi^2 xzi^^ ^) ' Pi^2 y{y)- it was discussed in Sec. [H instead of scanning the full Ar plane 
in this case it was sufficient to make two one-dimensional scans over Ax and Ay. In- 
stead of Ax, any other line in x-z plane may be taken, so for generality we consider the 
line inclined at an angle 6z from the x axis and denote the corresponding displacement as 
Ar;^ = (Ar^^ ■ cos6'^, 0, Ar^^ ■ sin6'^). 

If the spatial resolution also factorizes, V^^^{Sr) = V^f^{Sx, Sz) ■ Vy^^{6y), according to 
Eq. f|T5|l the vertex density can be split into two parts: 

7^ cose, d^K. ^^^^..^ ^ ^^^^^^^ ^^^^^ _ ^24) 



fATNiN2a dx dy dz 
the overlap integral from Eq. f|T5|) . its normalization follows from Eq. ffTOl) 

p^^(f, Ar) d^rd^Ar = 1. (25) 



Similarly, the relative normalization of and may be fixed by the requirement 

j pZf{x, z, Ar,,) dxdzdAr,, = j p^'^iy, Ay) dydAy = 1. (26) 
Following the same arguments as above one may then obtain the analogs of Eq. fpTl) : 

j pjf (x, z, Ar^^) dz2 dAr^^ ^ j '^^ ""^^^^ ~ ^X2)V^2{5x2) dSx2, 



pl^{y, Ay) dAy = j p^yiy - ^y)Vl^\6y) d6y, (27) 

where p2 x2 {p^i^)^ ^2 (Vy^^) are the second beam transverse profile along X2 (y) and the 
corresponding projection of the resolution, 

P2 X2{X2) = j Pf..(x, Z) dZ2, V,2{SX2) = J V^fiSx, Sz) d5z2 • (28) 

If one integrates Eq. 0241) over d^f, all information about the spatial distribution of 
vertexes is lost and the result should be expressible via rates. To be consistent with the 
normalization conditions f l26|l one should choose 



/" vx. A w ^ i^(Arg,, Ayo) 

/ p,.^ (x, z, Arz.^ dxdz = „^ , ^ — — - 

y ^ R{^r,,, Ayo)dl 



/i?(Ar^2, Ayo) dAr 

where Ar°^ and Ay^ are the fixed displacements during the scan performed along the other 
coordinate. As it was pointed out in Sec. [H they may be arbitrary. Indeed, since the rate 
is factorisable, Ar^^ or Ay^ appearing both in enumerator and in denominator cancel out. 
Note, that both Eqs. (127|) and (129|) lead to (l26l) after integration over X2, or over Ar^^^, 
Ayo, respectively. 

Taking into account Eqs. (1271) and (1291) . one can integrate (12^ over dydz2 dAr^z or over 
dx dz dAy and finally obtain 

• 6/iVvx(Ar^„Ayo)6/Ar^, i?(Arg„ Ayp) [ 

d^2 W = jRiArO^,Ay)dAy J ^^"^ " ^-^)^-(^-^) 

dNyj,iArl,Ay) dAy RjArl, Ayp) Z" M/1-Va, WA, 
Ty jRiAr..,Ay,)dAr.JP^y^'-''^''y ^'^^''^^ ^''^ 

where Nq = -^^J^^^^"^ = J Nyxi^f^z, ^y) dAr^zdAy. One can see that the distribution 
of vertexes accumulated during Ar^z i^y) scan, projected to X2 (y) axis, normalised and 
unfolded with the corresponding projection of the resolution Vx2 (Vy^^), gives the transverse 
beam profile p2 x2{x2) {p^i\{y))- 

If the scan is performed not continuously but stepwise, so that the data is taken at 
discrete points Ax = inej,.. Ay = ±mey with some integers n and m, the integrations may 
be approximated as discrete sums. This brings some systematic uncertainty, which may be 
estimated in the end from the reconstructed beam images. 
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3. Discussion and conclusions 

The collider luminosity and cross sections can be measured in van der Meer scan by 
sweeping the beams transversely across each other and by measuring the collision rate as 
a function of the beam displacement. We proved that the method remains applicable for 
the arbitrary beam crossing angle and the beam displacement plane, see Eqs. (|T0|) . (fT3|) and 

The four-dimensional beam densities may also be arbitrary. In particular, the formulas 
are valid for the beams initially mismatched in time. With non-zero crossing angle the time 
shift between the beams is equivalent to some shift in Ax, and therefore the integration 
over Ax in van der Meer method removes the dependence on timing. Also note, that in the 
case of x-y factorization, when it is sufficient to scan only along two perpendicular axes, 
their crossing point (Axq, Ayo) in Eq. f|T3|) is also arbitrary and not necessarily the point of 
maximal luminosity. The beams at this point can be mismatched in Ax or, equivalently, in 
time. 

The results are applicable to van der Meer scans performed at LHC accelerator in 2010. 
This was a primary tool for the absolute luminosity measurement at four major LHC exper- 
iments. The maximal crossing angle of 540 /xrad was in LHCb. It reduced the luminosity 
by about 4%, but caused a negligible correction, 7^ — 1 = 4 • 10~^, in van der Meer formula 

m- 

The original van der Meer method suggested in 1968 was based on counting the inter- 
actions. An important information is also contained in the spatial distribution of vertexes. 
With excellent modern detectors, like in LHC experiments, they are precisely measurable. 
We propose a new simple method how they can be used to reconstruct the beam images. 
The vertex distributions transverse to the beam and visible from its center should be accu- 
mulated during the scan and unfolded with the transverse vertex resolution. This should 
give the beam image in its transverse plane, see Eqs. (l2Ti) and (l30l) . The approach is valid 
for arbitrary beam shapes. 

From the reconstructed normalized beam profiles one can determine their overlap and 
then the luminosity using Eq. This gives an alternative way of the absolute luminosity 
measurement during van der Meer scan. Two methods are independent as the beam imaging 
method uses only the normalized shapes of the accumulated vertex distributions, while the 
traditional method uses only the integrals, i.e. the total number of interactions. 

There is one complication in determining the luminosity from the images, however. Ac- 
cording to Eq. ([2]), the luminosity depends on the overlap integral of the beam profiles in 
x-y laboratory plane, p^i^''^{r±)- With the proposed imaging method one can reconstruct 

the profiles pi2{^i2) transverse to the beams, see Eq. fl2T]) . They coincide with p^i2'^i^-i-) 
only for the coUinear beams, but in general lead to a different overlap integral. 

To correct for this effect, some information is needed on the distribution of particles 
along the beam. For example, in LHC one can do the following. As it was pointed out, 
the non-zero crossing angle reduced the luminosity and the overlap integral in 2010 scans 
by 4% or less. Since the effect is small, for its estimation it may be sufficient to approx- 
imate the bunch shapes along xi_2 and 2:1^2 as independent Gaussians with some effective 
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sigmas axi,2 and azi,2- In this case the x-projections in the laboratory frame are also Gaus- 
sian with sigmas cri2'^ = \/ (ca:i,2 cos 01,2)^ + (^"21,2 sin ai,2)^- If the bunch lengths of two 
beams are similar, a^i ~ crz2, they can be obtained from the width of the luminous re- 
gion (721,2 ~ V^'^ium- Here we assumed coUinear beams, but corrections due to non-zero 
crossing angle are negligible since in LHC cosai,2 ~ 1 and transverse sizes of the bunches 
are about three orders of magnitude smaller than their lengths. Since for Gaussian beams 
the luminosity is proportional to a/ crli + cr^2 5 its reduction due to the crossing angle is 
V<^li + ^x2/V(c^:ri cosai)2 + ((T^2 COS 02)^ + 2(af^jj^)2(sin^ «! + sin^ 02). 

If the detector is able to measure both the spatial coordinates and the time of the 
interactions, one can reconstruct not only transverse but the full three-dimensional beam 
images, see Eq. f l23p . The luminosity can be determined from them without any extra 
corrections or assumptions. 

The proposed imaging during van der Meer scan is very similar to the beam-gas imag- 
ing joj. The idea of the latter is to take a beam "photo" using the interactions with the 
gas remaining in the beam pipe, assuming it is distributed uniformly in the transverse 
plane. In both methods, after deconvolution with the vertex resolution, one can measure 
the transverse beam profiles and then the luminosity, taking into account the crossing angle 
correction. The beam-gas imaging method was successfully applied for the first time to 
measure the absolute luminosity in LHCb 01 . Its accuracy, as in van der Meer LHC scans in 
2010, was dominated by errors in bunch intensity Ni^2 measurement. The beam-gas method 
does not require moving of the beams and can be used during normal data taking. On the 
other hand, the advantage of the beam imaging during van der Meer scan is a much higher 
statistics of interactions and a localization of vertexes around a nominal luminous region 
where the vertex resolution is optimal. The methods have different systematics, and it is 
very advantageous to use both during van der Meer scan. 
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